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The effect of an external electric field F on the excitonic 
photoluminescence (PL) spectra of a symmetric coupled dou- 
ble quantum well (DQW) is investigated both theoretically 
and experimentally. We show that the variational method 
in a two-particle electron-hole wave function approximation 
gives a good agreement with measurements of PL on a narrow 
DQW in a wide interval of F including flat-band regime. The 
experimental data are presented for an MBE-grown DQW 
consisting of two 5 nm wide GaAs wells, separated by a 4 
monolayers (MLs) wide pure AlAs central barrier, and sand- 
wiched between Gao.7Alo.3As layers. The bias voltage is ap- 
plied along the growth direction. Spatially direct and indirect 
excitonic transitions are identified, and the radius of the ex- 
citon and squeezing of the exciton in the growth direction 
are evaluated variationally. The excitonic binding energies, 
recombination energies, oscillator strengths, and relative in- 
tensities of the transitions as functions of the applied field 
are calculated. Our analysis demonstrates that this simple 
model is applicable in case of narrow DQWs not just for a 
qualitative description of the PL peak positions but also for 
the estimation of their individual shapes and intensities. 

PACS numbers: 73.2Q.Dx, 78.55.-m, 78.66.-w, 78.66.Fd 



I. INTRODUCTION 

Because of their technical importance and unique phys- 
ical properties, semiconductor quantum wells have been 
the subject of intensive research since their first fabrica- 
tion in 1973. A brief review of the main achievements 
and a number of representative references can be found, 
e.g., in books ||l| and |^. The tunnel-coupled quantum 
states in DQWs are very sensitive to both electric and 
magnetic fields, and changes induced by these fields have 
been intensively studied by PL, PL excitation, and pho- 
toconductivity |,|,|,i3,|. 

A symmetric coupled DQW consists of two identical 
quantum wells separated by only a thin barrier. The 
energy levels of the coupled QWs split owing to the in- 
terwell tunneling. In the flat-band condition {F = 0), 
the eigenfunctions of the DQW have defined symmetries. 
In this situation, only transitions between electron and 



hole states of the same symmetry are optically allowed. 
The ground-state wave function is symmetric, while that 
of the first excited state is antisymmetric. When an elec- 
tric field is applied to the DQW, the wave functions be- 
come predominantly localized in one well only. Transi- 
tions allowed in the flat-band case evolve into spatially 
indirect (interwell) transitions, which shift linearly in en- 
ergy as a function of the applied field. On the other 
hand, transitions forbidden in the flat-band case become 
allowed and represent spatially direct (intrawell) tran- 
sitions. Two factors determine the optical transition in- 
tensities: (a) the overlap integral of the electron and hole 
single-particle wave functions, and (b) the population of 
the energy subbands for a given temperature. In order to 
describe the optical transitions correctly one has to ac- 
count for the Coulomb interaction between electrons and 
holes. We use a two-particle wave function composed of 
the electron and hole single-particle wave functions mul- 
tiplied by a function of their relative positions to get an 
exciton binding energy. 

The purpose of this report is to demonstrate that our 
quite simple theory, applied to the PL results obtained 
from our specifically designed sample, provides excellent 
quantitative agreement with regard to both the energetic 
position of the PL peaks and their shapes and intensities 
as a function of applied electric fields. These results rep- 
resent clearly a further improvement compared to the 
previously reported results ||,0,^,|| • 

II. THEORY 

The calculation starts from the envelope-function ex- 
citonic Hamiltonian 



Hex — Hoe + Hoh + H2D + U + E, 
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where the respective electron and hole single-particle 
terms read [Q 
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v — {e, h} denotes an electron or hole, and Qi, is the 
respective charge. mi,(zi,) and Vu{zu) are the respec- 



tive electron or hole Zy dependent effective masses and 
confining potentials of the DQW structure. The electric 
field F is applied along the z axis parallel to the growth 
direction. Hid represents the kinetic energy in the xy 
plane. The electron-hole exciton interaction is included 
by means of the Coulomb term [/, and the GaAs bandgap 
energy Eg completes the total energy H^x- 

Eigenenergies and eigenfunctions of the single-particle 
ID Hamiltonian Hqv are found by using linear combi- 
nations of analytical functions sin(^) and cos(^) (wells, 
F = 0), Airy functions Ai(^) and Bi(^) (wells, F > 0), or 
exp(±^) (barriers), and by matching the wave function 
amplitudes and their derivatives divided by the effective 
masses at each interface. Representative results of such a 
calculation made for a DQW structure with Lm = 18 ML, 
Lb = 4 ML (« 5.09 nm, and 1.13 nm, respectively), and 
for F = 30 kV/cm are schematically plotted in Fig. |l|. 



fied absolutely (ze, ^/i), only the relative distance between 
them projected onto the xy plane (p) is relevant. We 
minimize the total energy [Eq. (||)] for each exciton sepa- 
rately and obtain four wave functions {xn, X127 X21, X22}- 

In our calculation, we neglect the heavy- and light-hole 
bands mixing and assume strictly parabolic dispersion 
relations. Also, the mixing of the exciton states by the 
Coulomb interaction is not taken into account because 
of its small effect on the narrow DQW system pi. The 
GaAs/AlGaAs/AlAs material parameters and their com- 
positional dependence is taken as in p],p|ppC[|. 

The XijS are used to evaluate the optical oscillator 
strength (see, e.g., Ref. B) and the electron- hole overlap 
integral Fij(O), which is given by 

00 00 

Fij(O) = Xij{^e=z,zh=z,p=0)dz = Nij (p'i{z)^pj{z)dz. 
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FIG. 1. Scheme of the DQW potential profile with an ap- 
plied external field (solid lines). The energy levels (dotted 
lines) and wave functions (dashed lines) of electrons and holes 
are shown as well. Bases of the wave functions are at the cor- 
responding energies and their amplitudes are plotted in ar- 
bitrary units. The plot was calculated for L^ = 5 nm and 
F = 30 kV/cm. 

The first two single-particle eigenfunctions (p'({zy), 
f'^i^iy) of Hov are used to construct a set of basis func- 
tions for the variational calculation of the exciton states 
as a single product of one electron and one hole wave 
function multiplied by a function of relative electron hole 
positions ||J|: 



Xtj = N,jipl{z^)ip'l{zh)eyi-p 



y/p'^ + aij{ze - ZhY 

iti'j 



where the normahzation factors A^. 



(3) 
and the variational 



parameters Rij, atj are, in general, different for each ex- 
citon. While the z coordinates of both carriers are sped- 



(4) 

The normalized exciton PL intensity spectrum is then 
calculated as 
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[3 ~ l/ksT, kB is the Boltzmann constant and T the 
temperature. The convolution form Lij(E) expresses the 
normalized line-shape of the i,jth transition, which is 
determined by a thermal distribution of the excitonic ki- 
netic energy and slightly diffused by a Lorenz function. 
The width A^ — 0.5 meV is used in every case. 



III. EXPERIMENTAL 

Our PL experiments were performed on a sample 
grown by MBE at a temperature of 600°C on a semi- 
insulating GaAs substrate oriented in the [001] direc- 
tion. The growth started with a 500 nm wide n-doped 
(Si, 2 X IQi^ cm^^) GaAs layer, followed by a 300 nm 
wide n-doped (Si, 1.4 x 10^® cm"^) GaAlAs layer. The 
Al content in the GaAlAs layers was always 0.3. After 
this the following sequence of GaAlAs layers was grown: 
500 nm intrinsic, 5 nm p-(5-doped (C, 3 x 10^^ cm~'^) layer, 
and 100 nm intrinsic. On top of this separating layer, 
a sequence of three symmetric DQWs with 4 MLs (« 
1.13 nm) wide pure AlAs central barriers in between were 
grown, employing growth interruptions of 20 s at each 
heterointerface. The well widths are 35 MLs, 26 MLs, 
and 18 MLs (w 10 nm, 7.5 nm, and 5 nm), respectively. 
The DQWs are separated by a 100 nm wide GaAlAs layer 



in each case. The growth then continued with another 
sequence of GaAlAs layers: 100 nm intrinsic, 5 nm n-S- 
doped (Si, 4 x 10^^ cm~'^), and 500 nm intrinsic. On the 
top a 300 nm wide p-doped (C, 1.4 x 10^* cm"^) GaAlAs 
layer and a 20 nm wide p-doped (C,2 x 10^^ cm~'^) GaAs 
cap layer were grown. 

The p-i-n configuration of the sample allows us to ap- 
ply a bias voltage Upn by means of selective Ohmic con- 
tacts to the p-doped layers on the top and to the n-doped 
layers at the bottom. Usually, the bands in any p-i-n 
structure are tilted and application of forward bias is re- 
quired to flatten them. We, therefore, used a special sam- 
ple design comprising (5-doped layers inside the intrinsic 
region of the structure in order to screen the built-in 
electric field. Thus, the flat-band regime is obtained al- 
most at Upn = V, avoiding high dark current present 
at higher forward bias. The measured devices of the size 
250 /im x250 fim were defined photolithographically and 
mesa-isolated. Detailed results will only be presented for 
the narrow DQW with L^ = 5 nm and if, — 1.13 nm, 
for which the application of our simple theoretical model 
is justified. 
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FIG. 2. PL spectra of the DQW with 5 nm wide wells as a 
function of various bias voltages. The applied voltages Upn are 
denoted beside each curve. The inset shows the normalized 
integral intensity N/N{Up„ — OV) of the presented PL curves. 



Figure g shows the PL spectra of the sample for var- 
ious bias voltages. The sample was cooled in a closed- 
cycle cryostat and excited by a Ti:sapphire laser pumped 
by an Ar"'"-ion laser. We used an excitation power of 
~ 100 mW/cm^ at a photon energy of 1722 meV (below 
the bandgap energy of Gao.7Alo.3As). The emitted lumi- 
nescence was analyzed by a monochromator with 0.6 m 
focal length and a 1200 grooves/mm grating, and de- 
tected by a cooled CCD camera. We applied reverse bias 
up to Upn — —8 V with the maximum current of only 
~ 10 /iA. An estimation of a dissipated power and a neg- 
ligible shift in position of the spectrum measured with 
Upn — +2 V proved that Ohmic heating of the sample 
for both reverse and forward biased junction is insignifi- 
cant. 



IV. DISCUSSION 

The sample was designed as a single (not multiple) 
DQW structure in order to minimize the well width fluc- 
tuations. For the same reason, the growth interruptions 
were employed, and a pure AlAs barrier instead of an 
Ala:Gai_2;As one was chosen to improve its homogeneity 
and to avoid barrier height fluctuations. The broaden- 
ing of the direct exciton peak of ~ 4 meV (see Fig. g) 
and even lower at 10 K (~ 2 meV, not shown) indicates 
a quite good quality of the sample. Furthermore, the 
exciton peaks are not split due to large area monolayer 
fluctuations of the well widths, a phenomenon which is 
frequently observed on samples grown with interruptions. 

The (^-doped layers shielding the DQWs within the p-n 
junction are effective in case of an excitation below the 
Gao.7Alo.3As bandgap energy, when the carriers are ex- 
cited exclusively in the wells. The flat band condition is 
found almost at zero bias voltage. On the other hand, 
when using a He-Ne laser excitation above bandgap en- 
ergy, the (5-doped layers are neutralized and become in- 
active. Consequently, a forward bias corresponding to 
the bandgap energy is needed in order to overcome the 
built-in field and to establish the flat-band regime. 

The experimental spectra plotted in Fig. were ana- 
lyzed to gain maximum data for the theoretical proce- 
dure. In the calculation, we kept the parameters L^ and 
Lfc constant and varied only T and F for optimization. 
The temperature T = 45 K of the sample was obtained 
by a comparison of the measured and calculated relative 
peak intensities of the spectra in the linear field range. 
The field F was deduced by means of tracing the en- 
ergy distance between the first two exciton transitions, 
Xii and X12. A very good agreement of the calculated 
and observed transition energies is represented by Fig. y. 
Notice that the optimized T provides an excellent corre- 
spondence in the peak positions. The absolute positions 
of all the peaks including the resonant splitting of the 
symmetric-antisymmetric states were obtained without 
any free parameter. The small difference in the inten- 



sities and positions of the theoretical hnes ^ 1 meV in 
comparison with the experiment can be attributed to (i) 
the approximations used in the theory and (ii) small de- 
viations from the intended structure of the sample. 
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FIG. 3. PL peak positions as a function of electric field F. 
The lines are calculated for the DQW depicted in Fig. hi at 
T = 45 K. The markers indicate the peak positions of the 
experimental PL spectra plotted in Fig. H. 

The inset of Fig. ^ shows the deduced field F as a func- 
tion of the bias voltage applied to the sample. In the case 
of the spectrum close to flat-band regime (C/p„ = OF), the 
corresponding field [F = 1.6 kV/cm) was determined 
from a linear part of this dependence. As anticipated, 
the dependence is linear within a limited voltage inter- 
val {—iV < Upn < 21/), and saturates for higher reverse 
bias. The linear regime of the band tilting is in very good 
agreement with F calculated directly from the sample de- 
sign (dotted line in the inset of Fig. Q) . A small discrep- 
ancy can be explained as an effect of a non-homogeneous 
distribution of the electric field over the three DQW sys- 
tems, or as a result of a too simple description of the 
single wave functions in the applied variational method. 

Figure shows the spectra calculated for the studied 
DQW structure for the optimized temperature and the 
fields corresponding to individual PL curves plotted in 
Fig. 0. The spectra are magnified in accordance with 
the measured integral intensities N(Upn) (see the inset 
of Fig. S) to allow a comparison with the measured PL. 
The optimized temperature seems to be reasonable tak- 
ing into account that the sample in a closed-cycle cryo- 
stat is cooled by a cold finger and placed on a rather low 
thermally conductive substrate. The close concert of the 
experimental and calculated PL spectra in the low field 
interval (— 3^ < Upn < 2V) is evident. By increasing 
reverse bias the observed indirect transition diminishes. 
This results from a low stability of the lowest energy xn 
exciton which becomes dissociated by electron tunneling 
outside the DQW. The highest transition X22 appears 
near the flat-band case only. It manifests itself in Fig. 
in the spectrum calculated for F = 1.6 kV/cm as the 



small peak at the energy of 1626 meV close to the X21 
transition. In the experiment, these transitions merge 
into one peak observed at 1623 meV. Out of the flat-band 
case, X22 is suppressed due to a weak population of the 
higher levels and the lower value of the electron-hole over- 
lap integral i^22(0) [Eq. (Q)]. With respect to the good 
agreement of relative intensities of the spatially direct 
and indirect transitions we can conclude that the overlap 
integrals are well estimated and the particle distribution 
is described by our wave functions in a satisfactory way. 
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FIG. 4. Calculated PL spectra for the DQW depicted in 
Fig. hi at a temperature T = 45 K. The applied electric field 
F is given beside each curve. The inset shows the electric 
field induced in the DQW by the bias voltage Upn- 



V. CONCLUSION 

In this paper we have concentrated on the photolu- 
minescence of narrow symmetric coupled double quan- 
tum wells, where the energy difference between excitonic 
states is high enough to allow neglecting of band mix- 
ing effects. The inclusion of exciton interaction is un- 
avoidable for a correct description of the electron-hole 
states and positioning of optical transitions. We have 
shown that the simple model describes nearly perfectly 
all dominant features of the optical transitions, namely 



relative intensities, absolute positions of the individual 
transitions, and the resonant splitting of the symmetric- 
antisymmetric states. 
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